Attosecond streaking of correlated two-electron transitions in helium 
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We present fully ab initio simulations of attosecond streaking for ionization of helium accompanied 
by shake-up of the second electron. This process represents a prototypical case for strongly correlated 
electron dynamics on the attosecond timescale. We show that streaking spectroscopy can provide de- 
tailed information on the Eisenbud-Wigner-Smith time delay as well as on the infrared field dressing of 
both bound and continuum states. We find a novel contribution to the streaking delay that stems from 
the interplay of electron-electron and infrared-field interactions in the exit channel. We quantify all the 
contributions with attosecond precision and provide a benchmark for future experiments. 
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With the development of light-wave synthesis with at- 
tosecond precision time-resolved investigations of ultrafast 
electronic dynamics in atoms, molecules, and solids on 
the atomic time scale came into reach [1—3]. Attosecond 
streaking is one of the most fundamental processes in at- 
tosecond science allowing for a mapping of time informa- 
tion onto the energy axis yielding a time resolution in the 
order of a few (tens) of attoseconds [4-7]. It is a variant 
of a pump-probe setting with an ultrashort extreme ultravi- 
olet (XUV) pulse serving as pump and a phase-controlled 
few-cycle infrared (IR) field as probe [3], First proof-of- 
principle studies addressed the direct time-domain mea- 
surement of the life time of the Xe(4p _1 ) hole by Auger de- 
cay of ~ 8 fs [5] and the time -resolved photoemission from 
a tungsten surface by an energetic XUV pulse of ~ 300 as 
duration [6] providing insight into the relative time delay 
of ~ 110 as [6] (later ~ 85 as [8]) between photoemis- 
sion from core levels relative to the conduction band. Even 
shorter delays were more recently determined in the ele- 
mentary photoelectric effect where a time shift of the 2p 
relative to the 2s electron in neon as small as 21 as has been 
obtained in an attosecond streaking experiment by Schultze 
et al. [7]. These observations have triggered a flurry of 
theoretical investigations [7, 9-14]. While time delays of 
atomic photoionization on the one-electron (or mean-field) 
level as well as additional effective time shifts due to the 
dressing of the outgoing photoelectron [11, 12, 14] (of- 
ten referred to as "Coulomb-laser coupling" (CLC), [15]) 
could account for a delay of the order of 6 — 8 as, even the 
approximate inclusion of correlation effects failed to repro- 
duce the experimentally observed delay [7, 9]. However, 
multi-electron effects were only studied in stationary scat- 
tering approaches, i.e., by analyzing the phase of the dipole 
transition matrix elements for photoionization in the ab- 
sence of the streaking field. This group delay or Eisenbud- 
Wigner-Smith (EWS) delay (see e.g. [16]) is not directly 
what is measured by a streaking experiment [11]. It was 
suggested that the probing IR field might be responsible 
for the larger delay observed in the experiment [9]. There 



have been extensive time-dependent studies that simulated 
actual streaking experiments, but up to now only within the 
single-active electron approximation [11, 12]. The role of 
electron correlations and its interplay with the laser field 
has remained an open problem. Treating both the dressing 
of the atomic dynamics by the IR streaking field and the 
dynamical electron correlation exactly is still out of reach 
for neon but is possible for helium. 

In this letter we present streaking simulations for sin- 
gle ionization in helium including shake-up by solving 
the full two-electron time-dependent Schrodinger equation 
(TDSE) in the presence of the laser field without further 
approximation. In contrast to neon [7], where ionization 
from different atomic subshells was probed, shake-up ion- 
ization in helium relies strongly on electron-electron inter- 
actions. Here, the absorption of a single XUV photon leads 
to both emission of one electron and excitation of the other 
electron, which can only happen if the electrons share the 
photon energy and thus interact strongly. Studying ioniza- 
tion with and without shake-up may serve as benchmark 
for the interplay between mean-field, correlation, and field 
dressing effects. 

We show that for ionization without shake-up the streak- 
ing time shifts can be very well accounted for within the 
framework of Coulomb-laser coupling and the Eisenbud- 
Wigner-Smith delay. For ionization accompanied with ex- 
citation of the second electron correlation becomes impor- 
tant as field dressing effects modify the electron-electron 
interaction in the exit channel and give rise to an additional 
apparent time shift. We quantify all the contributions with 
attosecond precision and provide a benchmark for future 
experiments. 

In our computational approach (see e.g. [17, 18] for a 
detailed description) we solve the Schrodinger equation by 
the time-dependent close-coupling method, [19-21]. For 
the radial discretization we use a finite-element discrete- 
variable representation (FEDVR) [22-24], and propagate 
in time using the short-iterative Lanczos (SIL) algorithm 
[25, 26] with automatic time-stepping and error control. 
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FTG. 1. Single-ionization of He: (a) energy levels of (nl, el) final 
states; (b) electron energy distribution P~£(e, # = 0±10°) created 
by a 200 as (FWHM) XUV pulse with mean energy huj = 100 eV 
along the polarization direction with an opening angle of 10°. 
(c) streaking spectrogram of (b) with a 3 fs streaking pulse with 
A = 800 nm and intensity J ra = 4 • 10 11 W/cm 2 , lower (high mo- 



mentum) feature: 
up. 



shake-down, upper (low momentum): shake- 



The laser fields are linearly polarized and treated in dipole 
approximation. The XUV pulse has a Gaussian enve- 
lope with a FWHM duration of 200 as and an intensity of 
J XDV = 10 12 W/cm 2 (for which multiphoton processes can 
be neglected), while the 800 nm IR field has a sine-squared 
envelope with a FWHM duration of 3 fs and an intensity 
J IR = 4 • 10 11 W/cm 2 . For these laser parameters, numer- 
ically converged streaking simulations are obtained with a 
partial wave expansion with total angular momenta up to 
L = 8 and one-electron angular momenta up to /< = 5 and 
/> = 8. We use an asymmetric radial box with an extension 
up to 960 a. u. in one direction and 96 a. u. for the other 
radial coordinate (which is enough for ionic bound states 
up to n = 8 to be well represented). Each FEDVR element 
spans a length of 4.0 a.u. and contains a DVR of order 11. 
Atomic units are used throughout the paper. 

Fig. 1 shows the energy levels of helium and a corre- 
sponding single ionization spectrum including shake-up 
peaks ("correlation satellites"). Due to the large differ- 
ence in the binding energy between the He + (n = 1) and 
He + (n = 2) states (Fig. la) of 1.5 a.u. (40.8 eV) the two 



peaks are well separated and resolvable. For a typical XUV 
pulse (Huj = 100 eV) with 200 as duration (~ 9 eV spec- 
tral width), different shake-up channels (e.g. n = 2 and 
n = 3) are not resolved. However, numerically we have 
access to all channels separately. The electronic spectra 
for different delay times r of the ionizing XUV pulse rel- 
ative to the probing IR field build up a streaking spectro- 
gram (see Fig. lc), where the spectra are shifted in mo- 
mentum relative to the unperturbed asymptotic momentum 
p . The streaking spectrograms are obtained by projecting 
the propagated wavefunction ^{fx,fi,t) onto the single 
continuum constructed as a symmetrized product state of a 
bound state ; m (r) of the He + ion and a Coulomb wave 
V>k(r) with charge Z = l. 

The absolute time shifts i s are extracted by a nonlin- 
ear least-squares fit of the modified final momentum Pf(r) 
of the different channels (taken as the first moment of the 
electron spectrum) [11], 



p f (r) p - aA m (r + t s ) 



(1) 



where a is a correction factor for the amplitude of the mo- 
mentum shift induced by the streaking field. The result- 
ing t s contains information on the Eisenbud-Wigner-Smith 
(EWS) time delay £ews [27-29] of the atomic photoion- 
ization process in the absence of the streaking field as well 
as on the dressing of atomic and ionic states in the IR field. 
On the one-electron level the streaking time shifts ts can 
be decomposed as [11] 
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EWS 
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(2) 



The EWS time delay £ews is given by the energy deriva- 
tive of the phase of the dipole transition element (ipf \ z \ (j>o ) , 

t E ws(£) = ^arg ((^(£,0 = 0)1^)) , (3 ) 

which is evaluated along the laser polarization axis in for- 
ward direction, as is the streaking spectrogram. The second 
term £clc is an apparent time shift due to the interaction of 
the outgoing electron with both the long-ranged Coulomb 
potential and the IR streaking field [11, 12, 14] which ap- 
proximately scales with the final continuum energy E as 
~ —E~ 3/2 [11, 30]. We take t CLC (E) from the reference 
streaking shift ^ ls \E) of the pure Z = 1 Coulomb po- 
tential at the asymptotic electron energy E, 

t CLC (E) = tf ls \E) n 



dE 



o-i 



(4) 



where we have subtracted the EWS delay of the Coulomb 
phase <Ti in the final p continuum state for one-photon ion- 
ization of HQs). The CLC time shifts are only weakly de- 
pendent on the streaking laser field parameters, in particu- 
lar the wavelength of the IR field [12] and the duration of 
the XUV pulse. 

For strongly polarizable initial bound states another ap- 
parent time shift tisLC was found due to initial-state laser 
coupling [11, 31]. Energy shifts of the initial state due to 
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FIG. 2. Temporal shifts t s (■ : Is, □ : 2s, : 2p, • : n = 2) 
extracted from quantum mechanical streaking simulations and 
shifts predicted by Eq. 2, tcLC + ^EWS (lines, with £tslc = 
for He(ls 2 )) for single ionization of helium into an opening an- 
gle of 10° with respect to the polarization axis with and without 
shake-up of the second electron as a function of the final electron 
energy. Time shifts that belong to the same XUV energy are thus 
shifted by l[ 2 ^ - If' = 40.8 eV. Note that in the spectral re- 
gion of resonances (35 eV < E < 54.4 eV for an ionic Is state), 
streaking time shifts are not well defined. 

the interaction with the laser field are at the moment of 
ionization transferred to the final energy, and thus appear 
as apparent time shifts. 

A novel scenario appears for true multi-electron pro- 
cesses beyond the single-active electron (SAE) or mean- 
field level. We consider the prototypical two-electron pro- 
cess, the photoionization of He with shake-up/down, where 
electron correlation is expected to play a significant role. 
Fig. 2 shows the streaking time shifts i s for ionization with 
and without shake-up in comparison with the prediction 
*clc + *ews (Eq. 2). Initial state distortions tisLC are neg- 
ligible for the helium ground state He (Is 2 ). For shake- 
down to the ground state in the ionic system He + (ls), 
i.e., direct ionization without shake-up, the experimentally 
accessible streaking shift ts agrees remarkably well with 
Eq. 2. To calculate exact dipole matrix elements for sin- 
gle ionization with the correct boundary conditions even 
above the double ionization threshold (Eq. 3), we use the 
extraction method of Palacios et al. [32] based on exterior 
complex scaling and apply it to the wave packet z\cp ). For 
comparison we have also calculated the EWS time shift 
within the Hartree (mean-field) approximation where we 
take the ionic electron distribution to create an effective 
one-electron potential. When the bound electron is left 
in the ionic ground state, also the mean-field values t^ws 
agree with the exact £ews to within less than one attosec- 
ond (not shown). 

A surprisingly different picture emerges for shake-up to 
n = 2 (2s and 2p) where ts an d tcLC + ^ews strongly 
disagree (Fig. 2). The streaking time shifts for shake-up 



predict a delay with respect to shake-down for all inves- 
tigated photon energies which can not be accounted for 
by the corresponding EWS delays t^ws- Obviously, the 
interplay between electron-electron interaction and the IR 
streaking field strongly influences the streaking delay of 
the outgoing wavepacket. The additional time delay can 
be viewed as the result of the coupling of the dipole mo- 
ment of the shaken-up ionic state to the streaking field. 
Due to the entanglement of the system, the outgoing elec- 
tron still contains information on the ionic state. For the 
almost-degenerate n manifolds, this effect should be max- 
imal for intershell eigenstates of the dipole operator, i.e., 
Stark states (nk) with n = 2 and k = ±1. Indeed, the 
resulting time delays (see Fig. 3) are strongly enhanced 
(k = 1) and diminished (k = — 1) compared to ionic 
final states with well-defined angular momentum quan- 
tum number. There are two effects worth noting: First, 
the potential seen by the receding electron is modified by 
Vd{t) = d-f/r 3 , where d is the dipole moment of the 
Stark state. However, this time shift is already included 
in the (field-free) EWS delay and thus gives no additional 
measurement-induced time shift. Second, creation of a 
shaken-up electron in a Stark state in the presence of the IR 
field entails an additional energy contribution — d ■ F m (t) 
of the ionic system (see Fig. 3). This time-dependent en- 
ergy shift of the bound electronic final state of the resid- 
ual ion becomes thus visible in the energy and momentum 
distribution of the streaked ionized electron. Since this en- 
ergy shift is proportional to F m (t) it appears as an addi- 
tional time shift ij^c = a t an ( — ^mdf /k z ) /cjir. This 
process can be viewed as a two-electron generalization of 
the streaking shift for one-electron states with permanent 
dipole moment first discussed in [10]. For final ionic Stark 
states the streaking delay t s is indeed given to remark- 
able accuracy by the generalization of Eq. 2, 

ts = ^EWS + ^CLC + ^ISLC + ^FSLC (5) 

(with iisLc = for the He(ls 2 ) initial state). This 
streaking-field induced two-electron effect causes also the 
additional time shifts for final ionic states with well de- 
fined angular momentum, s and p, i.e., states with a neg- 
ligible dipole moment due to inversion symmetry. If we 
expand the final ionic states into Stark states, |2s/p) = 
l/y/2 (| 2+) ± 1 2— )) we can find an effective dipole mo- 
ment for the final two-electron states, 

< eff = (d + |c + | 2 + d-|c_| 2 )/(2| Ci | 2 ), (6) 

where c a = {2a, E6 \z\ipi). Note that the non-zero dipole 
moment in Eq. 6 is consistent with the inversion symmetry 
of the ionic bound state as this symmetry is broken by the 
selection of the emission direction of the streaked contin- 
uum electron. Even if the final state of the remaining ion 
is not detected, i.e., the time shift for the total n = 2 shell 
is measured, an averaged effect in the order of 5 as is visi- 
ble (see Fig. 3). Thus probing the momentum spectrum (or 
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FIG. 3. (a) Illustration of shake-up ionization with the ion re- 
maining in a Stark state and (b) the ionic energy in the laser 
field which leads to an additional momentum shift of the out- 
going electron, (c) Streaking time shifts tg (points) along the 
forward laser polarization direction for the two ionic parabolic 
states (n = 2, k = ±1) in comparison with n = 2, 2s and 2p states 
(see Fig. 2). In addition we show the complete prediction of is 
according to Eq. 5, icLC + *EWS + *fslc (solid lines). 



energy shift, respectively) of the ionized electron contains 
information on the dynamics of the remaining ion. 

We have thus identified an additional time shift iFSLC re- 
sulting from the back action of the excited bound state onto 
the continuum wave packet in the presence of the IR streak- 
ing field. Similar to the corresponding initial-state distor- 
tions, this additional time delay is not (or only weakly) de- 
pendent on the intensity of the IR field I m (within the range 
of useful intensities 10 10 W/cm 2 < J IR < 10 12 W/cm 2 ), 
cf. [11]. Thus, for photoionization of two-electron systems 
with shake-up, the total streaking time shift (Eq. 5) con- 
tains an additional contribution due to simultaneous exci- 
tation of the second electron. Attosecond streaking phases 
can therefore reveal information on dynamical polariza- 
tion in two-electron transitions. This information on the 
time-resolved bound-state excitation complements alterna- 
tive techniques of transient x-ray absorption spectroscopy 
[33] and transient bound state excitation by shaped two- 
color laser fields [34]. 

We conclude by discussing possible experimental ob- 
servations. While the streaking spectrogram for the final 
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FIG. 4. Relative total temporal shifts t^ n>1) 



(w) - t^\u) be- 



tween ionization with and without shake-up of the second elec- 
tron extracted from quantum mechanical streaking simulations 
for helium for different XUV energies, for shake-up into n = 2 
(blue line, open points) and sum of all shake-up states (n > 2, red 
line, solid points), lines to guide the eye. 



ground state, He + (ls), is well separated from that of two- 
electron excitation-ionization (Fig. 1), unambiguous obser- 
vation of the latter requires the separation of the n = 2 from 
higher shells with n > 3. This would require an XUV 
pulse with a Fourier-limited width corresponding to a du- 
ration of T X uv > 500 as (we have checked that our nu- 
merical results change within less than 1 as as compared to 
the 200 as used above). Alternatively, measuring Ly a pho- 
tons resulting from the radiative decay of the excited n = 2 
states in coincidence would allow to separate excitations of 
the n = 2 shell. Moreover, excitations of the 2s and 2p 
states could be separated by prompt vs. delayed Ly Q co- 
incidences in a weak DC electric field. A simple proof of 
principle experiment would be to monitor the large delay 
of shake-up into the n = 2 subshells, or an ensemble of ex- 
cited manifolds (n > 2), with respect to streaking without 
shake-up (n = 1) as a function of the XUV energy (Fig. 4). 
This would give rise to time advances up to 12 as when the 
n = 2 signal can be separated from higher shake-up states 
and up to 27 as when the total electron spectrum for n > 1 
is analyzed. The large difference in the delay for n = 2 
and n > 2 shows how sensitive the extracted streaking time 
shifts are to small contributions of different shake-up states 
in the streaking spectrum (see inset Fig. 1). The strong pho- 
ton energy dependence of the relative delay between n = 1 
and n > 1 is due to the variation of icLC with the final-state 
energy E. The ensemble of excited manifolds with n > 2 
would correspond to the streaking delay between the two 
peaks of the total ionization spectrum (see Fig. lb, black 
line) which should be relatively straight-forward and could 
serve as benchmark for the precision of the experimental 
streaking techniques. 

In summary, we have shown that for photoionization 
of helium streaking time shifts provide detailed informa- 
tion on the Eisenbud-Wigner-Smith time delay as long 
as Coulomb-laser coupling as well as laser-induced state 
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distortion effects are accounted for. For two-electron 
excitation-ionization the interplay of electron-electron and 
IR-field interaction in the exit channel leads to additional 
and novel contributions to the time shift of the outgoing 
wavepacket. We show that dynamical correlation effects 
can play a significant role for attosecond streaking exper- 
iments. Our theoretical ab-initio results can serve as an 
accurate benchmark for experimental attosecond streaking 
setups. 
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